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Abstract 

The saturation properties of neutron-rich matter are investigated in a relativistic mean-field 
formalism using two accurately calibrated models: NL3 and FSUGold. The saturation properties — 
density, binding energy per nucleon, and incompressibility coefficient — are calculated as a function 
of the neutron- proton asymmetry a = {N — Z)/A to all orders in a. Good agreement (at the 10% 
level or better) is found between these numerical calculations and analytic expansions that are 
given in terms of a handful of bulk parameters determined at saturation density. Using insights 
developed from the analytic approach and a general expression for the incompressibility coefficient 
of infinite neutron-rich matter, i.e., Kq^o) = Kq + Kroi^ -|- . . ., we construct a Hybrid model with 
values for Kq and Kj^ as suggested by recent experimental findings. Whereas the Hybrid model 
provides a better description of the measured distribution of isoscalar monopole strength in the 
Sn-isotopes relative to both NL3 and FSUGold, it significantly underestimates the distribution of 
strength in ^'^^Pb. Thus, we conclude that the incompressibility coefficient of neutron-rich matter 
remains an important open problem. 
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I. INTRODUCTION 



The incompressibility of neutron-rich matter remains at the forefront of both experimental 
and theoretical investigations due to of its fundamental role in constraining the equation 
of state (EOS) of cold dense matter. The incompressibility coefficient characterizes the 
small-density fluctuations of infinite nuclear matter around its equilibrium point, thereby 
providing the first glimpse into the "stiffness" of the equation of state. 

It is widely accepted that the Giant Monopole Resonance (GMR) — the quintessential 
compressional mode — provides the cleanest, most direct route to the nuclear incompress- 
ibility. In a procedure first proposed by Blaizot and collaborators IH |2], finite nuclei GMR 
energies as well as the nuclear-matter incompressibility should both be computed within 
the same theoretical framework. This procedure avoids altogether the reliance on macro- 
scopic ( "liquid-drop-like") approaches that have proven unreliable for the extraction of the 
incompressibility coefficient of infinite nuclear matter from finite nuclei GMR energies [21 S] ■ 
In most theoretical approaches, an accurately calibrated model is obtained by fitting the 
model parameters to a set of selected ground-state properties of finite nuclei. In some recent 
instances, excited states — computed as the consistent linear response of the mean-field 
ground state — have also been incorporated into the fit |5]. If such accurately calibrated 
models are able to reproduce the experimental distribution of monopole strength (or at 
least some of its moments) then the value of the incompressibility coefficient predicted by 
the model is regarded as reliable. Following this procedure it has been established that the 
incompressibility coefficient of symmetric nuclear matter falls within the following relatively 
narrow range: 7^0 = 240 ±10 MeV PE]. 

Because of the collective nature of the GMR, a strong coherent peak develops only in 
the case of relatively heavy nuclei. Indeed, the monopole strength in "light" nuclei, such 
as ^°Ca, is strongly fragmented. As (stable) heavy nuclei are characterized by a significant 
neutron excess, experimental studies of the GMR probe the incompressibility of neutron-ricli 
matter rather than that of symmetric matter. As such, GMR energies on a variety of nuclei 
having different neutron-proton asymmetries [a = {N — Z)/A], such as ^^Zr, ^^^Sn, ^^*^Sm, 
and ^°*Pb, provide simultaneous constraints on the incompressibility of symmetric nuclear 
matter (Kq) as well as on its leading a correction, a quantity that will be denoted by Kt-. 
That is, the incompressibility of infinite neutron-rich matter may be parametrized to leading 
order in the neutron-proton asymmetry as Ko{a) = Kq + K-rO^. It is therefore natural to 
assume that previous lessons learned in the case of Kq will remain relevant for K-r. First 
and foremost, should not be inferred from an extrapolation to the oo limit from 
laboratory experiments on finite nuclei. Rather, one should continue to follow the procedure 
advocated by Blaizot and demand that the values of both Kq and K^^ be those predicted by 
a consistent theoretical model that successfully reproduces the experimental GMR energies 
of a variety of nuclei. We note that in the present contribution both Kq and K^- refer to 
bulk properties of the infinite system. 

It is therefore the aim of the present manuscript to: (a) use accurately calibrated relativis- 
tic mean-field models to extract the saturation properties of infinite neutron-rich matter; 
(b) compute GMR energies for a variety of nuclei using the consistent isoscalar-monopole 
response of the mean-field ground state; and (c) confront these theoretical results against 
experimental GMR energies — especially the new data on Tin [HI [7]. As a byproduct of this 
procedure, analytic approaches to the saturation properties of infinite neutron-rich matter 
based on a few bulk parameters calculated at saturation density will be validated against 
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exact numerical results. 

The manuscript has been organized as follows. In Sec. |Tl]a brief review of the relativistic 
mean-field formalism will be provided. Particular emphasis will be placed on developing a 
thorough description of the properties of infinite neutron-rich matter and on the various bulk 



parameters that define its behavior around nuclear-matter saturation density. In Sec. Ill 
results will be presented for the evolution of the saturation point with neutron-proton asym- 
metry using both exact (numerical) and approximate (analytic) approaches. We finish this 
section by revisiting the topic of why is Tin so soft? [6l [3 [HI [9l [10]. Our summary and 
conclusions will follow in Sec. [IVl 



II. FORMALISM 

The starting point for the calculation of various ground-state properties is an interacting 
Lagrangian density of the following form: 



^int = ^ [gs4>- [g^V^+^T ■ b^+|(l+r3)A^) 7^ 



- ^W)^-^(^7s0)^ + ^(<?^V;V^) Va.(^?>, ■ b^) [glV.V^) . (1) 

The Lagrangian density includes an isodoublet nucleon field {ip) interacting via the exchange 
of two isoscalar mesons, a scalar (0) and a vector iV^), one isovector meson {h^), and the 
photon {A^) [Til Il2]- In addition to meson-nucleon interactions, the Lagrangian density is 
supplemented by nonlinear meson interactions with coupling constants denoted by k. A, C, 
and Ay that are responsible for a softening of the nuclear-matter equation of state, both 
for symmetric-nuclear and pure-neutron matter. Whereas k and A are instrumental in the 
softening of the equation of state of symmetric nuclear matter near saturation density, C, 
softens the equation of state but at higher densities. Finally, the mixed vector coupling 
Av [13] has been introduced to soften the density dependence of the symmetry energy — a 
quantity that is predicted to be stiff in most relativistic mean-field models. This effective 
Lagrangian has been used to compute a variety of ground-state observables at the mean-field 
level and will be used here to study the incompressibility of neutron-rich matter. Further 
details on the calibration and implementation of the relativistic mean-field models may be 
found in Refs. [H] |T2l [ISl [E] and references therein. 

Asymmetric nuclear matter is an idealized system consisting of an infinite number of 
neutrons and protons interacting exclusively through the nuclear force. At zero temperature 
and in the thermodynamic limit (where both the baryon number A = N + Z and the 
volume of the system V tend to infinity) the binding energy per nucleon depends solely 
on two intensive variables, the baryon density p = A/V and the neutron- proton asymmetry 
a = {N — Z)/A (the latter may be expressed as a = {pn — Pp)/p in terms of the nucleon 
densities). By studying such an idealized system, one hopes to elucidate how the volume 
and symmetry terms of the semi-empirical mass formula \TE\ evolve with density. 

It has become customary to write the energy per particle of infinite nuclear matter as 
follows: 

E/Aip, a)-M = Sip, a) = ^snm(p) + a'52(p) + a^S^{p) + . . . , (2) 

where we have indicated that E/A is measured relative to the nucleon rest mass M. As 
the neutron-proton asymmetry is constrained to the interval < a < 1, the total energy 
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per particle E{p.a) is often expanded in a power series in a^. Note that odd powers of 
a do not contribute to the expansion owing to the symmetry of the strong force between 
hke-nucleon pairs. The leading term in this expansion, i.e., ^snm(p) =^{P: o: = 0), represents 
the contribution from symmetric {N — Z = A/2) matter. The leading C(Q;^)-correction to 
the symmetric limit, i.e., S2{p) =S{p), is known as the symmetry energy. The contribution 
a'^S{p) thus measures the energy involved in converting part of the protons in symmetric nu- 
clear matter to excess neutrons, at total baryon density p. The above power-series expansion 
is particularly useful as the symmetry energy dominates the corrections to the symmetric 
limit for all values of a. Indeed, to an excellent approximation the energy per particle of 
pure neutron matter {a = l) may be written as follows: 

^^PNM (p) = S{p,a = l) !V £sNM (p) +S{p) . (3) 

The main feature that we aim to understand in the present manuscript is the evolution 
with neutron-proton asymmetry of the bulk properties of infinite nuclear matter — such as 
the saturation density, the binding energy at saturation, and the incompressibility coeffi- 
cient. Particularly important is to characterize the sensitivity of the results to the density 
dependence of the symmetry energy. To do so and to establish a baseline, we start by 
describing the behavior of symmetric nuclear matter near saturation density. 



A. Symmetric Nuclear Matter 

One of the hallmarks of the nuclear dynamics is the saturation of symmetric (a = 0) 
infinite nuclear matter. The saturation point is characterized by an equilibrium density of 
about p„ ^0.15 fm^"^ and an energy per particle of e„ ~ — 16 MeV. Given that the pressure 
P = p^ dSsNu/dp of symmetric nuclear matter vanishes at saturation, then the small density 
fluctuations around the saturation point are fully characterized by the incompressibility 
coefficient Kq. Yet for reasons that will become clear later, the behavior of symmetric nuclear 
matter is expanded in a Taylor series up to third order in the small density fluctuations. 
That is, 

1 1 
^^snm(p) = £o + 2^ox^ + Q^ox^ + ■ ■ ■ ' 

where x is a conveniently defined dimensionless parameter that characterizes the deviations 
of the density from its saturation value, i.e., 

x-^. (5) 

The various bulk coefficients that characterize the behavior of the symmetric system near 
saturation density are given as follows: 

^0 = ^snm{x = 0) = i^sNM(p = Po) ) (6a) 



SNM 



dp'^ 
dp^ 



(6b) 

P=Po 



Qo = 4V(^ = 0) = 27p^('%M) . (6c) 



P=Po 
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B. Symmetry Energy 



In the so-called "parabolic" approximation the deviations from the symmetric {a = 0) 
limit are controlled by the (9(a^)-symmetry energy [see Eq. As has been done for 

the symmetric case, the behavior of the symmetry energy around nuclear-matter saturation 
density may be conveniently characterized in terms of a few bulk parameters, namely, 

S{p) = J + LX+ ^i^symX^ + ^QsymX^ + . . . , (7) 

where J, L, Kgy^, and Qsym are the values of the symmetry energy, slope, curvature, and 
third derivative at saturation density. However, unlike symmetric nuclear matter the "sym- 
metry pressure" L does not vanish. Indeed, the symmetry pressure — a quantity that 
strongly influences the neutron-skin thickness of heavy nuclei — is (within the parabolic 
approximation) directly proportional to the pressure of pure neutron matter. That is, 

^0 = IpoL . (8) 

A one-parameter fit to the low- density behavior of the symmetry energy that is frequently 
used in the heavy-ion community is of the following form [161 [13 HI [13 [20l |2T1 [22] : 

\ 7 



S{p)=So(^] =J{l + 3x)\ (9) 

where in arriving to the last term we have made use of Eqs. ^ and ([T]). To the extent that 
the above parametrization is accurate, the following relations should be satisfied: 

^^(iL^^^- (lOa) 
Ksym = (1^) = 9 J7(7 - 1) , (10b) 

Po = PoJl- (10c) 

With due caution, mainly because the connection of heavy-ion experiments to the EOS 
is not at all trivial and often involves model-dependent extrapolations of the measured 
data, significant constraints on the value of the coefficient 7 have been extracted in the 
last few years from different experimental observables. For instance, in intermediate-energy 
heavy-ion collisions, the analysis of isoscaling data [201 EI] provides a 7 value around 0.69. 
Transport model simulations of data related to isospin diffusion favor the milder constraint 
7 ~ 0.69-1.05 [m [22]. Some nuclear collective modes provide another tool to probe the 
behavior of S{p) at subsaturation densities. The values Pq = 2.3 ± 0.8 MeV/fm^ and 
J = 32 ± 1.8 MeV extracted in Ref. [23] from pygmy dipole resonances suggest a value 
of 7 ~ 0.5 ± 0.15, whereas the constraint 23.3 < S{p = 0.1 fm"^) < 24.9 MeV obtained 
in Ref. [21] from the properties of the giant dipole resonance in ^°^Pb hints at a value of 
7 ~ 0.5-0.65. These findings from experimental isospin-sensitive signals imply a rather soft 
nuclear symmetry energy at subsaturation densities. An analysis [25] of a set of neutron 
skins of nuclei measured across the mass table by antiprotonic techniques yields a similar 
conclusion. Finally, recent studies of the low-density behavior of pure neutron matter using 
universal properties of dilute Fermi gases seem to support the same findings [26l [2Z] ■ 
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C. Neutron- Rich Matter 



Insofar as neutron-rich matter saturates, the energy per particle may continue to be 
written as in Eq. Q. Thus, the aim of this section is to characterize the evolution of 
the saturation point — in particular, the saturation density, the binding energy per nucleon, 
and the incompressibility coefficient — as a function of the neutron-proton asymmetry a = 
(N — Z)/A. To do so, it will prove instructive to proceed along two alternative paths: 
one purely analytical and the other purely numerical. In the analytic case, the saturation 
properties of neutron-rich matter will be derived from a handful of bulk parameters that 
characterize the behavior of both symmetric nuclear matter and the symmetry energy around 
saturation density, as was done in Eqs. ^ and ([T]). This purely analytic procedure, already 
well known in the literature, will be contrasted against a numerical procedure that is free 
of any assumptions or approximations beyond that of the mean-field approximation. We 
will verify that these two alternative approaches agree at the few percent level, thereby 
lending support to the analytic approach in elucidating the evolution of the incompressibility 
coefficient with neutron excess. 

According to Eqs. ^ and ([T]), the energy per particle of asymmetric nuclear matter with 
a neutron-proton asymmetry a may be written in the form 

S{p, a) ^ {e, + Ja^) + La^x + ^{Kq + a^K^ym)x^ + ^(<5o + a^Qsyrn)x'^ ■ (11) 

Clearly, the presence of the linear term shifts the saturation point from xo = to xq, where 
the latter is defined as the physical solution to the following equation: 

— \= a^L + {Ko + a^K,yra)x + - (Qo + a'Qsym)a;' = . (12) 

Although the roots of this equation may be found by solving a simple quadratic equation, 
the 0(q;^) solution may be solved by inspection. One obtains 

x, = -^a^ + 0{a^) or p /p^ = 1 + = 1 - 3^a^ + C(a^) . (13) 

The values for the energy per particle and the incompressibility coefficient may now be found 



by expanding Eq. (11) around this new (a^O) value of the saturation density: 

S{p,a) = S{x^,a) + ^{x - x,,fS"{x^,a) + . . . . (14) 
Alternatively, by introducing the dimensionless parameter 

- ^ (15) 

to characterize the deviations of the density from the new saturation point, the expansion 
of £{p,a) given in Eq. (14) may be cast in a form analogous to Eq. ^ for £^snm(p)- That 
is, 
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£^(p, a) = £^(Xo,a) + - (1 + 3a;J £^"(xQ,a) x + . . . = e^{a) + -KQ{a)x +..., (16) 
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Model 


ms 






9^ 


9'. 






A 


c 


Av 


FSU 


491.500 


782.500 


763.000 


112.1996 


204.5469 


138.4701 


1.4203 


+0.02376 


0.06 


0.03 


NL3 


508.194 


782.501 


763.000 


104.3871 


165.5854 


79.6000 


3.8599 


-0.01591 


0.00 


0.00 


Hybrid 


508.194 


782.501 


763.000 


106.2575 


165.5848 


79.6483 


4.5472 


-0.01952 


0.00 


0.00 



TABLE I: Model parameters used in the calculations. The parameter k and the meson masses ms, 
rriv, and nip are all given in MeV. The value of the nucleon mass is taken as M = 939 MeV. 



where the energy per particle and the incompressibihty coefficient at the new saturation 
density are given by 



£{xg, a) = Eg + Ja^ + 0{a^ 



(17a) 
(17b) 



The analytic results correct to second order in a are summarized in the following set 
of equations where the quantities pr, Er, and Kr represent the deviations of the saturation 
density, energy per particle, and incompressibihty coefficient of infinite matter away from 
the symmetric N = Z limit: 



Po(") = Po + Prtt^ + C'(a'^) = Po - 3Po^"^ + , 
eg{a) = Eg + ErO^ + 0{a^) = e,, + Ja^ + 0{a^) , 

Ko{a) = Ko + K^a" + 0{a^) = + (i^sym - 6L - + O(a') 



(18a) 
:i8b) 
(18c) 



In view of the profuse choices of terminology existing in the literature, our notation and 
conventions are discussed further in the Appendix. 

III. RESULTS 



Having developed the necessary formalism in the preceding sections, we devote this sec- 
tion to present the results of our calculations. As we have done elsewhere |H], our results 
were generated using two accurately calibrated models: NL3 [281 EH] and FSUGold [5]. In 
addition, we performed calculations with a Hybrid model to be introduced later. Effective 
meson masses (i.e., interaction ranges) and coupling constants for the present models are 
displayed in Table [l] as defined by the Lagrangian density of Eq. ([T| . 

With the above sets of parameters, one may compute the nuclear-matter equation of state, 
namely, the energy per particle as a function of density and neutron excess. In particular, 
one can extract values for the various bulk parameters defined in Eqs. ^ and ([t]) that 
characterize the behavior of neutron-rich matter around saturation density; these parameters 
are listed in Table |TTj Note that in order to make contact with the parametrization given 
in Eq. ([9]), the value of the exponent 7 listed in Table [H] was extracted from a fit to the 
symmetry energy in the low-density range of p = (0.3 — 1.0)pp. The found results are in 
consonance with the prediction 7 = L/3J that follows from Eq. ^ (namely, 7 = 0.62 for 
FSUGold and 7 = 1.06 in the case of NL3 and the Hybrid model). Whereas the FSUGold 
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Model 


Po 


eo 




Qo 


J 


L 


K 


7 


FSU 


0.148 


-16.30 


230.0 


-523.4 


32.59 


60.5 


-51.3 


0.64 


NL3 


0.148 


-16.24 


271.5 


+204.2 


37.29 


118.2 


+100.9 


0.98 


Hybrid 


0.148 


-16.24 


230.0 


-71.5 


37.30 


118.6 


+110.9 


0.98 



TABLE II: Bulk parameters (as described in the text) characterizing the energy of symmetric 
nuclear matter [Eq. Q] and the symmetry energy [Eq. ([7|] at saturation density. All quantities 
are in MeV, with the exception of po given in fm~^ and the dimensionless parameter 7 defined in 
Eq. ^. 



and NL3 models agree on the energy and density at saturation — quantities that are tightly 
constrained by existent ground-state observables — significant discrepancies emerge in all 
remaining parameters. The main difference between the two models may be succinctly 
summarized by stating that whereas FSUGold predicts a soft behavior for both symmetric 
nuclear matter (through Kq) and the symmetry energy (through L), NL3 predicts a stiff 
behavior for both. Note that "stiff" or "soft" refers to whether the energy increases rapidly 
or slowly with density. 

That symmetric nuclear matter and the symmetry energy are either both soft (as in the 
FSUGold model) or both stiff (as in the NL3 model) may lie at the core of the problem 
in reproducing the experimentally measured OMR energies in the Sn-isotopes [HI El |8]. 
According to Eq. (18c), a large value of L (as in NL3) generates a large softening of the 
incompressibility coefficient relative to its value in symmetric nuclear matter. However, a 
large incompressibility coefficient in symmetric nuclear matter Kq (as in NL3) hinders the 
softening generated by K^. Conversely, FSUGold predicts a relatively small value for Kq. 
However, its soft symmetry energy generates a small (absolute) value for Kj. that precludes 
the significant reduction in the incompressibility coefficient required by the experimental 
GMR energies. In an effort to circumvent this problem — and this problem only — we have 
generated a Hybrid model having the same incompressibility coefficient as FSUGold while 
preserving the stiff symmetry energy of NL3 (see Table |ll]). As seen in Table |l( this was 
accomplished through a slight adjustment of the scalar self-coupling parameters k and A. 
Note that it is not our intent to accurately calibrate the Hybrid model introduced here. 
Rather, the Hybrid model — although reasonably accurate — should merely be regarded 
as a "test" model that illustrates how surprisingly soft are the experimental GMR energies 
of the Tin isotopes relative to the theoretical predictions. Indeed, as we will display later 
in Fig. [8} not even such an artificially-tuned model can fully account for the rapid softening 
of the GMR energies with neutron excess. Let us mention that nowadays it is generally 
acknowledged that experimental data on compressional-mode giant resonances support a 
value of Kq ~ 240 MeV [30] for the incompressibility coefficient of symmetric nuclear matter. 
Recently measured data on the breathing mode of Sn isotopes seem to favor a constraint 

= —550 ±100 MeV for the asymmetry term in the nuclear incompressibility [SI [7]. A 
similar value of ~ —500 MeV has been reported from independent experimental evidence 
available from isospin diffusion observables in nuclear reactions [IHl [22] and from neutron 
skins of nuclei [25] . Tables [H] and HI confirm that the Hybrid model is consistent with both 
of the indicated Kq and K^ values. 

The evolution of the equation of state with neutron-proton asymmetry a is displayed in 
Fig. [ijfor the three models considered in the text: FSUGold (solid blue lines), NL3 (dashed 
green lines), and Hybrid (dot-dashed red lines). The a = curve corresponds to symmetric 
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' p/p„ ^ 

FIG. 1: (Color online) Density dependence of the energy per particle in infinite nuclear matter 
at the shown neutron-proton asymmetries according to the FSUGold (blue solid line), NL3 (green 
dashed line), and Hybrid (red dot-dashed line) models. 



nuclear matter whereas the a = l curve corresponds to pure neutron matter. In all models 
nuclear matter ceases to saturate at a value of the neutron-proton asymmetry slightly larger 
than a = 0.75. Figure |2] provides an expanded version of the results for symmetric nuclear 
matter [a = 0) and for pure neutron matter {a = 1). To an excellent approximation — 
especially at the subsaturation densities of relevance to this work — the difference between the 
equation of state of pure neutron matter and that of symmetric matter equals the symmetry 
energy of Fig. 3j On the other hand, a truncated expansion S{p) = J + Lx + ^Ksymx'^ 
[cf. Eq. ([t])] appears to provide a fair enough representation of the actual value of S{p) in 
a range of densities roughly between half and twice the saturation density of symmetric 
nuclear matter. Indeed, we find that the discrepancies of this approximation compared with 
the exact S{p) are less than 5% in a density range 0.45p(, < p < 2.7p,) for FSUGold and 
0-33po < p ^ 2.15Pq for NL3 and the Hybrid model. 

One observes that in all cases the Hybrid model seems to follow closely the predictions of 
the NL3 model in Figs. [l}]3| For the symmetry energy this has been done by construction. 
For symmetric nuclear matter, however, the Hybrid model is indeed softer than NL3 — and 
as soft as FSUGold — at saturation density. That the Hybrid model tracks NL3 at high 
density is a reflection of the vector self-coupling parameter ( having been set to zero in 
both models. This conflrms that the value of the incompressibility coefficient of symmetric 
nuclear matter at saturation density has practically no impact on the equation of state of 
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high-density matter and, by extension, on most neutron-star properties [31|. 




From the evolution of the equation of state with a we conclude quite generally that as the 
neutron-proton asymmetry increases, the saturation density is reduced, the binding energy 
weakens, and the nuclear incompressibility softens. Based on the structure of Eqs. (18), we 
regard these trends as model independent. Indeed, that the saturation density is reduced 
[i.e., p,-<0) is a simple reflection of (a) symmetric nuclear matter has a stable equilibrium 
point [Ki^ > 0) and (b) the pressure of pure neutron matter at saturation density is positive 
(L > 0). Further, that the binding energy weakens {i. e., er> 0) follows from the fact that pure 
neutron matter is not self bound, namely, J> |eo|. Finally, that the nuclear incompressibility 
softens requires Kj. < 0. This is the hardest condition to satisfy as it depends on higher- 
order derivatives of the equation of state, namely, on Qo and -K'sym [see Eq. (18c)]. However, 
barring anomalously large values for these two quantities, the condition Kr < hinges also 
on the pressure of pure neutron matter at saturation density being positive. This is due to 
the large coefficient multiplying L in Eq. (18c), which provides the dominant contribution 
to Kr as compared to Ksjm and Qq/Kq. We conclude that whereas the signs of p,-, ^t, 
and Kr are fairly model independent, their model-dependent magnitude is determined by 
two fundamental parameters of the equation of state: the incompressibility coefficient of 
symmetric nuclear matter -ft'o>0 and the symmetry pressure L. 

To examine the evolution of the saturation point with neutron-proton asymmetry we 
have tabulated in Table III the values for p,-, Er-, and Kr- These quantities (which have 
been enclosed in parenthesis) were computed directly from the analytic expressions given in 
Eqs. (18). The bulk parameters that they depend on were previously extracted from a fit to 
the equation of state of symmetric nuclear matter and to the symmetry energy and are listed 
in Table ITT] As alluded earlier, the fact that neither Qq nor -ft'sym are anomalously large in 
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FIG. 3: (Color online) Density dependence of the symmetry energy coefficient for the investigated 
models. 



Model 




er(MeV) 


i^^(MeV) 


FSU 
NL3 
Hybrid 


-0.117 (-0.117) 
-0.188 (-0.194) 
-0.215 (-0.229) 


32.60 (32.59) 
37.24 (37.29) 
37.20 (37.30) 


-275.45 (-276.77) 
-682.65 (-697.36) 
-531.98 (-563.86) 



TABLE III: Leading 0(a^)-correction to the evolution of the saturation density, energy per particle, 
and incompressibility coefficient of asymmetric nuclear matter. Quantities outside the parentheses 
were extracted from a quadratic fit to the numerical results in the 0<a'^<0.1 range, whereas the 



quantities in parenthesis were computed from the analytic expressions given in Eqs. (18) 



the present models, results in an asymmetry term in the nuclear incompressibility that 
is dominated by the symmetry pressure L. In particular, note that the value of in the 
Hybrid model is consistent with the value extracted in Refs. [6l [7] from the measurement of 
the GMR energies in the Sn-isotopes. 

But how accurate are the expressions given in Eqs. (18)? To test the reliability of these 
analytic expressions we have carried out a purely numerical exercise that is exact within the 
purview of the mean-field approximation. There is no reliance on the parabolic approxima- 
tion, as in Eq. ([2]), or on expansions around the saturation density of symmetric matter, such 
as in Eq. (11). Basically, the equation of state of asymmetric nuclear matter is computed 
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numerically for a range of values of the neutron-proton asymmetry in the < < 0.1 range. 
For each value of a, the new saturation point — namely the density, energy-per-nucleon, and 
incompressibility coefficient at the minimum — is computed. Once this procedure is com- 
pleted, one extracts the three desired coefficients (p,-, e^, and K^-) from a least-squares fit 
to the numerical data. Such a procedure is illustrated in Figs. |4| [5} and |6] for the saturation 
density, energy-per-particle, and incompressibility coefficient, respectively. In all cases the 
inset includes a comparison between the exact numerical results (displayed with lines) and 
the analytic approximations (displayed with symbols). Moreover, pr, Er, and Kr have also 
been tabulated in Table |III[ The agreement between the analytic and numerical results is 
fairly good (at worst the discrepancies amount to ~ 6%) suggesting that arguments based on 
Eqs. (18), which imply an expansion in both a and the density around the saturation point 
oi N = Z matter, are not only insightful but also quantitatively accurate. The goodness 
of the parabolic approximation for the binding energy of asymmetric matter seems to be 
confirmed also in other frameworks such as microscopic many-body calculations [321 EH] or 
model analyses of the symmetry energy coefficient in nucleus-nucleus collisions 
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FIG. 4: (Color online) Evolution with increasing neutron-proton asymmetry of the baryon density 
that corresponds to the exact saturation point of asymmetric nuclear matter. The inset displays by 
symbols the result of a least-squares fit in the 0<a2<0.1 range assuming a parabolic dependence. 

We finish this section by revisiting a topic recently addressed in the literature: why is 
Tin so soft? [SI [71 El m [10]. Namely, GMR energies of even- A isotopes of Tin from A = 112 
to A = 124 measured in a recent experiment [6l [7] are significantly lower than the values 
predicted with accurately calibrated, otherwise successful, mean field models. Note that the 
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same models satisfactorily predict the GMR excitation energy of ^*^Zr and ^°^Pb. Therefore 
we pose the following question: can the Hybrid model succeed where the other two (FSUGold 
and NL3) have failed? Recall that the Hybrid model was built with the explicit purpose of 
having a "low" incompressibility coefficient of Kq ^ 230 MeV and a "large" (and negative) 
asymmetric term of K^- = —532 MeV (see Table HI), unlike FSUGold where both Kq and 
\Kr\ are low, and unlike NL3 where both Ko and \Kr\ are high. Thus, as in Ref. [8], the 
distribution of isoscalar monopole strength for the even-A Tin isotopes — from ^^^Sn up to 
^^■^Sn — was computed in a relativistic random-phase approximation (RPA). Details of the 
method may be found in Ref. [35] . 
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FIG. 5: (Color online) Same as described in the caption to Fig. |4]for the energy per particle. 



Figures [7], [8| and [9] are reminiscent of those published in Ref. [S| but in the present case 
results are also included for the Hybrid model. It is clear from Fig. [7] that the experimental 
distribution of strength in the Tin isotopes is best reproduced by the Hybrid model. Note 
that RPA distributions of strength fail to account for the full — escape-plus-spreading — 
width of the resonance. Whereas the RPA calculation accounts properly for the escape width 
{i.e., the coupling to the continuum is treated exactly) it fails to account for its spreading 
component as this one is related to configurations significantly more complex than those 
included in the RPA approach. 

On the other hand, the RPA approach is sophisticated enough to reproduce the experi- 
mental centroid energy of the resonance. As it was done experimentally, the centroid energy 
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FIG. 6: (Color online) Same as described in the caption to Fig. [4] for the incompressibility coeffi- 
cient. 



was computed from the ratio of the mi to the mo moment. That is, 



where S\,[(i^lS) is the distribution of strength. The integration Umits have been fixed at 
LO\ = 10 MeV and UJ2 = 20 MeV, respectively, and the integrals have been evaluated at the 
small momentum transfer of g~0.23 fm~^ (or g~45 MeV) [8]. The theoretical predictions 
for -Eqmr in the Tin isotopes are displayed in Fig. [8] in comparison with the experimental 
data from RCNP ^61 17| and TAMU ^301 [Ml EZ]- 

Although the FSUGold and Hybrid models share the same value of the incompressibility 
coefficient in symmetric nuclear matter, the Hybrid model provides a softer distribution of 



strength because of its largest (negative) asymmetric term Kr (see Table HI). Ultimately, 
this result hinges on the fact that the Hybrid model has, as NL3, a significantly stiffer 
symmetry energy. All in all, the agreement between the Hybrid model and experiment is 
adequate, although the model — indeed all models — could benefit from a steeper slope in 
the change of the centroid energy Eqmr with mass number A. To test the robustness of our 
results we have used an improved version of the Hybrid model that was obtained through 
a slight adjustment of the scalar mass m^ (508.194 494 MeV) and the corresponding 
coupling constant g"^ (106.2575 100.4048, which yields the same gl/ml value and ensures 
that all of the properties of the EOS of infinite nuclear matter remain unaltered). This 
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FIG. 7: (Color online) Comparison between the distribution of isoscalar monopole strength in all 
neutron-even ^^^Sn-^^^Sn isotopes extracted from experiment j6l [7j (black solid squares) and the 
theoretical predictions of the FSUGold (blue solid line) , NL3 (green dashed line) , and Hybrid (red 
dot-dashed line) models. 
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mild adjustment yields better ground-state masses for a few selected nuclei (^°Ca, ^°Zr, and 
^°^Pb), albeit at the expense of slightly worsen charge radii. Yet the GMR energies for the 
Sn-isotopes get softened by at most 1.5%. This confirms one of the main results of this 
work, namely, that even a model with a soft Kq (such as FSUGold) and a stiff Kt- (such as 
NL3) is unable to fully account for the rapid softening of the experimental GMR energies 
in the Tin isotopes. 
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FIG. 8: (Color online) Comparison between the GMR centroid energies (mi/mo) in all neutron- 
even ^^^Sn-^-^^Sn isotopes extracted from experiment |Sl[Tj (black solid squares) and the theoretical 
predictions of the FSUGold (blue up-triangles), NL3 (green down-triangles), and Hybrid (red dot- 
dashed line) models. Also shown (filled gold circles) are experimental results from the Texas A&M 
group [301 EH [37] for the cases of "^Sn, ^^^Sn, and ^^^Sn. 



Therefore, where does theory stand with respect to experiment? Unquestionably, a Hy- 
brid model having a soft incompressibility coefficient but a stiff symmetry energy leads to a 
significant improvement when compared with the experiment on the Tin isotopes [HI E] . Un- 
fortunately, the Hybrid model does not fare as well against other observables (see Table IV ) . 
First, the Hybrid model predicts a GMR centroid energy in ^osp]-, i?QMR = 13.27 MeV (or 
13.16 MeV if we use the Hybrid model with rris = 494 MeV), significantly lower than the 
experimental value of 14.17±0.28 MeV [30]; in contrast, the FSUGold model gives a value 
of £^GMR = 14.04 MeV [5] that is consistent with experiment. Note that similar trends have 
recently been reported by Avdeenkov and collaborators [10]. To appreciate the significant 
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FIG. 9: (Color online) Distribution of isoscalar monopole strength for ^°Zr, ^^^Sn, ^^^Sm, and 
^'^'^Pb as predicted for the FSUGold (thick lines) and Hybrid (thin lines) models. Experimental 
centroid energies for these nuclei were reported in Ref. |30] and have been tabulated along the 



theoretical predictions in Table IV 



softening of the Hybrid model relative to FSUGold, the distribution of isoscalar monopole 
strength for ^°Zr, ^^^Sn, ^^^Sm, and ^°'^Pb is displayed in Fig. [o] Second, as we have argued 
earlier, a large negative asymmetry term Kr requires a large positive value of the asym- 
metry pressure L. However, models with a stiff symmetry energy appear to be in conflict 
with model-independent predictions for the equation of state of pure neutron matter at low 
densities [2S1 123 ESI EH] ■ Finally, note that a stiff symmetry energy at densities below sat- 
uration also seems to be disfavored by heavy- ion data [IS], [211 [22], although in this case 
the model dependence may be more significant and the results are not always without some 
controversy jlQl [H] . 



IV. CONCLUSIONS 

The saturation properties of neutron-rich matter were studied as a function of the 
neutron-proton asymmetry within the framework of relativistic mean-field models. We ob- 
served that infinite nuclear matter continues to saturate up to values of the neutron-proton 
asymmetry of the order of a = (A^ — Z)/A< 0.75. Moreover, it was found quite generally 
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- 26 
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14.17 ±0.28 


14.04 


14.32 


13.27 



TABLE IV: Giant-Monopole-Resonance centroid energies (Sgmr = '^i/"T'o in MeV) obtained from 
the distribution of monopole strength integrated over the ijJ\-lo2 range for those nuclei studied in 
Ref. ^30j; a represents their respective neutron-proton asymmetry. 



that as infinite nuclear matter departs from the symmetric (A^ = Z) limit, the saturation 
density lowers, the binding energy weakens, and the nuclear incompressibility softens. 

The manuscript was organized around three main themes: (a) the use of accurately 
calibrated relativistic mean-field models to extract the saturation properties of neutron-rich 
matter directly from numerical computations and the comparison of these numerical results 
against approximate analytic approaches; (b) the use of the same models to compute the 
distribution of isoscalar-monopole strength in various nuclei; and (c) the comparison of these 
theoretical predictions against the experimentally measured GMR energies. 

To make contact between the equation of state of bulk neutron-rich matter and GMR en- 
ergies on finite nuclei, the incompressibility of infinite neutron-rich matter was parametrized 
in terms of two bulk parameters, namely, Kq and K^^ with the former being the incompress- 
ibility coefficient of symmetric matter and the latter parametrizing the (small) deviations 
from the symmetric limit [see Eq. (18c)]. Note that never in the manuscript we relied 
on semi-empirical {liquid-drop-like) formulas to extract properties of infinite-matter from 
extrapolating finite-nuclei results to the oo limit. In this manner we followed the time- 
honored tradition initiated by Blaizot and collaborators [T| 12] of demanding that the values 
of both Kq and Kr be those extracted from a consistent theoretical model that successfully 
reproduces the experimental GMR energies of a variety of nuclei. 

As part of the first theme, the evolution with neutron-proton asymmetry of the satu- 
ration density, binding energy per nucleon, and incompressibility coefficient were extracted 
from a fit to the numerically generated equation of state. Once these properties were ex- 
tracted, their dependence on the neutron-proton asymmetry a was captured through a 
simple parametrization in powers of with no reliance on the parabolic approximation 
of Eq. Q nor on an expansion involving bulk-model parameters, as in Eq. ((ill). Having 
completed this numerical procedure, we explored the possibility of reproducing the exact 
numerical results from analytic expansions based on a few bulk parameters of the equation 
of state determined at normal nuclear-matter saturation density [see Eqs. (18)]. For all 
three bulk properties — the saturation density, the binding energy per nucleon, and the 
incompressibility coefficient — the analytic values were in close agreement with those com- 
puted numerically. This seems to be a robust result as it holds for all three (FSUGold, NL3, 
and Hybrid) models; see Table HI Thus, we concluded that the analytic expressions are 
not only insightful but also quantitatively accurate. Particularly interesting was the case of 
the incompressibility coefficient K.^ that is given as the sum of three potentially ''large and 
cancelling" contributions. However, we found that one of these three terms — the slope of 
the symmetry energy L — dominates Kr [see Eq. (18c)], thereby making sensitive cancel- 
lations unlikely. This result revealed an interesting correlation between Kr and L that may 
be further explored by the upcoming Parity Radius Experiment (PREx) at the Jefferson 
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laboratory. PREx promises to measure the neutron radius of ^"^^Pb accurately and model 
independently via parity- violating electron scattering |l3l HI]. PREx will provide a unique 
experimental constraint on the density dependence of the symmetry energy due its strong 
correlation to the neutron radius (or neutron skin thickness) of heavy nuclei [U] . 

To test the predictions of these three models contact had to be made with available 
experimental data on finite nuclei. Thus, the distribution of isoscalar monopole strength 
was computed for a variety of nuclei in a consistent RPA approach [35j- In particular, the 
main motivation behind introducing the Hybrid model was due to the inability of both 
FSUGold and NL3 to reproduce the recently-measured GMR energies along the isotopic 
chain in Tin [6[ [7] . By adopting a relatively small value for the incompressibility coefficient 
in symmetric matter [Kq = 230 MeV) together with a fairly large negative value for its 
leading deviation from the symmetric limit [K^- ^ —530 MeV), we were constructing a Hybrid 
model with a significantly softer incompressibility coefficient for neutron-rich matter. Such 
a softening indeed produced a significant improvement vis-a-vis the experimental data on 
the Sn-isotopes; see Fig. [8j Whereas FSUGold and NL3 overestimate the centroid energy 
in ^^^Sn by about 0.7 and 1.0 MeV, respectively, the Hybrid model falls within 0.1 MeV of 
the experimental data. Indeed, the predictions of the Hybrid model fall within 0.1 MeV of 
the experimental data for the full isotopic chain if one takes account of the uncertainties in 
the data. However, although the improvement in the case of the Sn-isotopes is significant 
and unquestionable, an important problem remains: the Hybrid model underestimates the 
GMR centroid energy in ^''^Pb — the heaviest doubly-magic nucleus — by almost 1 MeV. 
This suggests that the rapid softening with neutron excess predicted by the Hybrid model 
may be unrealistic. 

Thus, where does theory stand with respect to experiment? One possibility, given that 
FSUGold reproduces the centroid energy in both ^°Zr (with a = 0.11) and ^osp]^ (with 
a = 0.21), is that its predictions for Kq and Kr are reliable but that its failure to reproduce 
the GMR energies in Tin is due to missing physics unrelated to the incompressibility of 
neutron-rich matter. We feel inclined to favor this possibility because of two main reasons. 
First, the missing physics may be to some extent related to the open-shell structure of 
the Tin-isotopes, a property that makes pairing correlations essential and endows Tin with 
its superfiuid character. Support in favor of this scenario has been recently presented in 
Ref. |16] where a surface pairing force was used to bring theory much closer to experiment, 
at least from ^^^Sn to ^^°Sn. Second, the large and negative value suggested from the 
experimental extraction of Kr may be at odds with theoretical constraints deduced from 
the behavior of dilute Fermi gases that seem to suggest a moderate value for the pressure 
of pure neutron matter at saturation density [27]. (Note that the pressure of pure neutron 
matter, or equivalently the slope of the symmetry energy L, largely determines the behavior 
of Kt-.) This suggests that the value of ii'T- = — 550±100 MeV inferred from experiment [HI Ej 
may suffer from the same ambiguities already encountered in earlier attempts to extract the 
incompressibility coefficient of infinite matter from finite-nuclei extrapolations. Yet the final 
resolution on "why is Tin so soft?" awaits further theoretical insights. 
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APPENDIX A: NOMENCLATURE AND TERMINOLOGY 

In this section we address what we perceive as a confusing state of affairs in regard to 
the nomenclature used to characterize the symmetry energy. First, we note that no uniform 
terminology exists even to denote the neutron-proton asymmetry coefficient. Indeed, the 
symbols / [47J , a [221 ill |19] , /3 |32] , and 5 IIS |50] are all used in the literature to denote 
the neutron-proton asymmetry coefficient [N — Z)/A of asymmetric nuclear matter. Second, 
and perhaps even more confusing, is the myriad of different symbols used to refer to the 
same bulk properties. For example, all of the following expressions may be found in the 
literature [iTlltH]: 

S{p) = J + Lx + ^-KsymX^ + ... (Ala) 



2 



2 



04 + ^(p - Po) + T^ip - Pof + • • • > (Ale) 



where x = {p—pg)/3pg. Moreover, another plausible expansion of the symmetry energy may 
be around the equilibrium Fermi momentum [50]. That is, 

S{p) = J + Ly+ ^K^y^y^ + ... , (A2) 

where the deviation from the equilibrium Fermi momentum has been parametrized in terms 
of the dimensionless parameter y defined as 

(A3) 

Recall that the Fermi momentum and the baryon density are related by the following ex- 
pression: 

One potential confusion between the two different expansions of the symmetry energy (in 
terms of either x or y) is that in the presence of a linear term (such as L) the coefficients 
are in general not equal. Indeed, the various bulk coefficients are related as follows: 

J = J , (A5a) 
L = L, (A5b) 

-^sym = -^sym + 2L . (A5c) 
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That is, at order and higher, the expansion coefficients in terms of the Fermi momentum 
y differ from the corresponding ones used in an expansion in terms of the density x. 

It is also common practice to express the finite nucleus incompressibihty coefficient (Ka) 
by means of a hquid-drop-hke mass formula [21 EJ HI SHI EOl EB E3, which highlights the 
physical meaning of the various contributions to Ky^. That is, 

Ka = i^vol + i^surf^-'/' + ( ) + ^CoulJ^ + ... . (A6) 

In some works, the coefficient Kt- is denoted by f^asy [IHl 1221 SH] or K^s [HI E2]- To add to 
the confusion in notation, in the original contributions by Blaizot and collaborators [H |21 EO] 
the term Kgym was used instead of K^- in Eq. (A6). It appears that at present -ft'sym has 



been "universally" adopted to refer to the curvature of the symmetry energy at saturation 
density, as in Eq. ([T]). 

In summary, we adopt the following convention in the present manuscript — and hopefully 
in all future works. The energy-per-particle of asymmetric nuclear matter is denoted as 
follows: 

S{p, a) = SsMp) + S{p)a'' + 0{a^) 

= (^e, + ^Kox^ + ^Qox^ + . . . ) 

+ (^J + LX+ ^i^symX' + ^QsymX^ + + 0{a^) , (A7) 

where the dimensionless parameters x and a characterize the deviations from saturation 
density and from the symmetric limit, respectively. That is, 

(A8) 

Finally, the quantities pr, Sr, and Kr, have been introduced to denote, respectively, the 
evolution with the neutron-proton asymmetry of the saturation density, the energy per 
particle, and the incompressibihty coefficient of infinite neutron-rich matter. That is, 

p,{a) = p, + pra^ + 0{a^) , (AlO) 
e, (a) =€,+ Era^ + C(a^) , (All) 
Ko{a) =Ko + K^a'^ + 0{a^) . (A12) 



21 



J. p. Blaizot, D. Gogny, and B. Grammaticos, Nucl. Phys. A265, 315 (1976). 

J. P. Blaizot, J. F. Berger, J. Dcchargc, and M. Girod, Nucl. Phys. A591, 435 (1995). 

J. M. Pearson, Phys. Lett. B271, 12 (1991). 

S. Shlomo and D. H. Youngblood, Phys. Rev. C47, 529 (1993). 

B. G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett 95, 122501 (2005), nucl-th/0504034. 
T. Li et al., Phys. Rev. Lett. 99, 162503 (2007), 0709.0567. 
U. Garg ct al., Nucl. Phys. A788, 36 (2007), nucl-ex/0608007. 
J. Piekarewicz, Phys. Rev. C76, 031301 (2007), 0705.1491. 

H. Sagawa, S. Yoshida, G.-M. Zeng, J.-Z. Gu, and X.-Z. Zhang, Phys. Rev. C76, 034327 
(2007), 0706.0966. 

A. Avdccnkov et al. (2008), 0808.0478. 

B. D. Scrot and J. D. Walccka, Adv. Nucl. Phys. 16, 1 (1986). 

B. D. Serot and J. D. Walecka, Int. J. Mod. Phys. E6, 515 (1997), nucl-th/9701058. 

C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 5647 (2001), astro-ph/0010227. 
B. G. Todd and J. Piekarewicz, Phys. Rev. C67, 044317 (2003), nucl-th/0301092. 

A. Bohr and B. R. Mottelson, Nuclear Structure (World Scientific Publishing Company, New 

Jersey, 1998), vol. I: Single-particle motion. 

M. B. Tsang et al., Phys. Rev. Lett. 86, 5023 (2001), nucl-ex/0103010. 
M. B. Tsang et al., Phys. Rev. Lett. 92, 062701 (2004). 

D. V. Shetty et al., Phys. Rev. C70, 011601 (2004), nucl-ex/0406008. 

B. -A. Li and L.-W. Chen, Phys. Rev. C72, 064611 (2005), nucl-th/0508024. 

D. V. Shetty, S. J. Yennello, and G. A. Souliotis, Phys. Rev. C75, 034602 (2007), nucl- 
ex/0505011. 

D. V. Shetty, S. J. Yennello, and G. A. Souliotis, Phys. Rev. C76, 024606 (2007), 0704.0471. 
B.-A. Li, L.-W. Chen, and C. M. Ko, Phys. Rept. 464, 113 (2008), 0804.3580. 
A. Klimkiewicz et al., Phys. Rev. C76, 051603 (2007). 

L. Trippa, G. Colo, and E. Vigezzi, Phys. Rev. C77, 061304 (2008), 0802.3658. 

M. Centelles, X. Roca-Maza, X. Vinas, and M. Warda, Phys. Rev. Lett. 102, 122502 (2009), 

0806.2886. 

A. Schwenk and C. J. Pethick, Phys. Rev. Lett. 95, 160401 (2005), nucl-th/0506042. 
J. Piekarewicz, Phys. Rev. C76, 064310 (2007), 0709.2699. 

G. A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C55, 540 (1997), nucl-th/9607039. 

G. A. Lalazissis, S. Raman, and P. Ring, At. Data Nucl. Data Tables 71, 1 (1999). 
D. H. Youngblood, H. L. Clark, and Y.-W. Lui, Phys. Rev. Lett. 82, 691 (1999). 

H. Mueller and B. D. Serot, Nucl. Phys. A606, 508 (1996), nucl-th/9603037. 

W. Zuo, I. Bombaci, and U. Lombardo, Phys. Rev. C60, 024605 (1999), nucl-th/0102035. 

H. Hcisclbcrg and M. Hjorth-Jensen, Phys. Rept. 328, 237 (2000), nucl-th/9902033. 

S. K. Samaddar, J. N. De, X. Vinas, and M. Centelles, Phys. Rev. C78, 034607 (2008), 

0804.2459. 

J. Piekarewicz, Phys. Rev. C64, 024307 (2001), nucl-th/0103016. 

D. H. Youngblood et al., Phys. Rev. C69, 034315 (2004). 

Y. W. Lui, D. H. Youngblood, Y. Tokimoto, H. L. Clark, and B. John, Phys. Rev. C70, 
014307 (2004). 



22 



C. J. Horowitz and A. Schwenk, Nucl. Phys. A776, 55 (2006), nucl-th/0507033. 
S. Kowalski et al., Phys. Rev. C75, 014601 (2007), nucl-ex/0602023. 

D. V. Shetty et al. (2006), nucl-ex/0603016. 

A. Ono, P. Danielewicz, W. A. Friedman, W. G. Lynch, and M. B. Tsang (2005), nucl- 
ex/0507018. 

B. K. Agrawal, S. Shlomo, and V. K. Au, Phys. Rev. C68, 031304 (2003), nucl-th/0308042. 

C. J. Horowitz, S. J. Pollock, P. A. Souder, and R. Michaels, Phys. Rev. C63, 025501 (2001), 
nucl-th/9912038. 



R. Michaels, P. A. Souder, and G. M. Urciuoh (2005), URL |http :77iiallaweb . j lab . org/ 
|parity/prex[ 

B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 

J. Li, G. Colo, and J. Meng, Phys. Rev. C78, 064304 (2008), 0811.4091. 

V. Baran, M. Colonna, V. Greco, and M. Di Toro, Phys. Rept. 410, 335 (2005), nucl- 

th/0412060. 

R. J. Furnstahl, Nucl. Phys. A706, 85 (2002), nucl-th/0112085. 

L.-W. Chen, C. M. Ko, and B.-A. Li, Phys. Rev. C76, 054316 (2007), 0709.0900. 

J. P. Blaizot, Phys. Rept. 64, 171 (1980). 

T. V. Chossy and W. Stocker, Phys. Rev. C 56, 2518 (1997). 

S. K. Patra, M. Centelles, X. Vinas, and M. Del Estal, Phys. Rev. C65, 044304 (2002), 
nucl-th/0111007. 



23 



